Series expansions of the density of states in SU{2) lattice gauge theory 
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We calculate numerically the density of states n{S) for SU{2) lattice gauge theory on L'' lattices. 
Small volume dependence are resolved for small values of S. We compare ln{n{S)) with weak 
and strong coupling expansions. Intermediate order expansions show a good overlap for values 
of S corresponding to the crossover. We relate the convergence of these expansions to those of 
the average plaquette. We show that when known logarithmic singularities are subtracted from 
ln{n{S)), expansions in Legendre polynomials appear to converge and could be suitable to determine 
the Fisher's zeros of the partition function. 
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I. INTRODUCTION 

Quantum Chromodynamics is a widely accepted the- 
ory of strong interactions. From a theoretical point of 
view, understanding the large distance behavior in terms 
of the weakly coupled short distance theory has been an 
important challenge. The connection between the two 
regimes can be addressed meaningfully using the lattice 
formulation. In the pure gauge theory (no quarks) de- 
scribed with the standard Wilson's action, no phase tran- 
sition between the weak and strong coupling regime has 
been found numerically for SU{2) or SU{3) and the the- 
ory should be in the confining phase for all values of 
the coupling. Recently, convincing arguments have been 
given pj, [2] in favor of the smoothness of the renormaliza- 
tion group flows between the two fixed points of interest, 
putting the confining picture on more solid mathematical 
ground. 

The absence of phase transition discussed above sug- 
gests that it is possible to match the weak coupling and 
the strong coupling expansions of the lattice formula- 
tion. However, if we consider these two expansions, for 
instance for the average SU{2) plaquette as a function of 
f3 = 4/5^, we see in Fig. [T]that there is a crossover region 
(approximately 1.5 < (3 < 2.5) where none of the two ex- 
pansions seem to work. This behavior is probably related 
to singularities in the complex (3 plane [H, |3| that are not 
completely understood. In the case of the one plaque- 
tte model [b\ , taking the inverse Laplace transform with 
respect to (3 (Borel transform) of the partition function 
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FIG. 1: Weak and strong coupling expansions of the aver- 
age plaquette P for SU{2) at various orders in the weak and 
strong coupling expansion compared to the numerical values. 



yields a function that has better convergence properties. 
It would be interesting to know if this feature persists on 
V = L'^ lattices. 

In this article, we study expansions of the inverse 
Laplace transform of the partition function (the density 
of states) of SU{2) lattice gauge theory on symmetric 
4 dimensional lattices. The density of states is denoted 
n{S) and defined precisely in Sec. |IT1 It gives a rela- 
tive measure of the number of ways to get a value S of 
the action. Knowing n{S), we can calculate the parti- 
tion function and its derivatives for any real or complex 
value of /3. In particular, it could be used to determine 
the Fisher's zeros of the partition function [1, 0, [!]■ The 
choice of SU (2) is motivated by the existence of a par- 
ticular symmetry Q which allows to determine the be- 
havior of n{S) near its maximal argument without extra 
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calculation. In Sec. IIIIl we explain why ln{n{S)) is ex- 
pected to scale like the volume and can be interpreted 
as a "color entropy". Numerical calculations of ri(5) ob- 
tained by patching plaquette distributions multiplied by 
the inverse Boltzmann weight at values of /3 increasing by 
a small increment are presented in Sec. IIVI The article 
is focused on comparisons with numerical data on a 6* 
lattice where finite volume effects are not too large and 
plaquette distributions broad enough to allow a smooth 
patching. The values of n{S) on such lattice are com- 
pared with those on a 4** and 8^ lattice. It is interesting 
to note that the volume dependence is resolvable only for 
small values of S where a behavior ln{S)/V is observed 
for ln{n{S)) 

The numerical results are compared with expansions 
that can be obtained from the strong (Sec. |V| and weak 
(Sec. I VI I) coupling expansions of the average plaque- 
tte. Intermediate orders in these expansions show a good 
overlap for values of S that correspond to the crossover. 
We then show that the convergence of the new series 
can be related empirically to those of the series for the 
average plaquette. The weak coupling expansion deter- 
mines the logarithmic singularities of ln{n{S)) at both 
boundaries. When these singularities are subtracted we 
obtain a bell-shaped function that can be approximated 
very well by Legendre polynomials (Sec. IVII[) . We con- 
clude with possible applications for the calculations of 
the Fisher's zeros and open problems. 



II. THE DENSITY OF STATES 

We consider the standard pure gauge partition func- 
tion 



^ = n / '^'^'^^ 



PS 



with the Wilson action 

5 = ^(1 - il/N)ReTr{Up)) 



(1) 



(2) 



and (3 = 2N/g'^. We use a D dimensional cubic lattice 
with periodic boundary conditions. For a symmetric lat- 
tice with sites, the number of plaquettes is 

Mp = L^D{D - l)/2 . (3) 

In the following, we restrict the discussion to the group 
SU{2) and D = 4. For SU{2), one can show 0] that 
the maximal value of S is 2Afp. We define the average 
plaquette: 



P = {S/Afp) = -d{lnZ/Up)/df3 . 



(4) 



Inserting 1 as the integral of delta function over the 
numerical values 5 of 5 in Z , we can write 



Z= dSn{S)e 
Jo 



with 

n{S) = n / '^'^'^(^ - - i'^/N)ReTr{Up))) (6) 

We call n{S) the density of states . A more general dis- 
cussion for spin models [T^j or gauge theories 11 1 can 
be found in the literature where the density of states is 
sometimes called the spectral density. From its defini- 
tion, it is clear that n{S) is positive. Assuming that the 
Haar measure for the links is normalized to 1, the par- 
tition function at /3 = is 1 and consequently we can 
normalize n{S) as a probability density. 

A first idea regarding the convergence properties of 
various expansions can be obtained from the single pla- 
quette model 0- In that case, we have 



lip; 



(S) = -^S{2-S) . 



(7) 



(3S 



(5) 



The large /3 behavior of the partition function is deter- 
mined by the behavior of n(S) near S — 0. In this ex- 
ample, n{S) cx ^/S for small S, implies that Z oc /J""^/^ 
at leading order. Successive subleading corrections can 
be calculated by expanding the remaining factor \/2 — S 
in powers of S and integrating over S from to cxd. If 
we factor out the leading behavior, we obtain a power 
series in 1/(3. The large order behavior of this power 
series is determined by the large order behavior of the 
expansion of \^2 — S, itself dictated by the branch cut 
at S* = 2. One can see Q that the S- integration over 
the whole positive real axis converts an expansion with a 
finite radius of convergence into one with a zero radius 
of convergence. On the other hand, if the S'-integration 
is carried over the interval [0, 2], the resulting series con- 
verges but the coefficients need to be expressed in terms 
of the incomplete gamma function. From this example, 
one may believe that it is easier to approximate n{S) than 
the corresponding partition function. However, it is not 
clear that these considerations will survive the infinite 
volume limit. Note also that the behavior of n{S) near 
S = 2 can be probed by taking /3 ^ — oo in agreement 
with the common wisdom that the large order behavior 
of weak coupling series can be understood in terms of the 
behavior at small negative coupling. 

It was showed @ that if the lattice has even number of 
sites in each direction and if the gauge group contains —1, 
that it is possible to change j3ReTrUp into —f3ReTrUp 
by a change of variables Ui ^ —Ui on a set of links such 
that for any plaquette, exactly one link of the set belongs 
to that plaquette. This implies 

Z{-P) = e'^'^^'-Zip) (8) 

This symmetry implies that 

n{2Mp -S)= n{S) . (9) 

In the following, we will be working exclusively with 
SU{2) which contains —1 and lattices with even num- 
bers of sites in every direction. We will thus assume that 
Eq. ^ is satisfied and we only need to know n{S) for 

< 5 < Tv;. 
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III. VOLUME DEPENDENCE 

In this section, we discuss the volume dependence of 
the density of state. We make this dependence exphcit 
by writing n{S,Afp). Given the density of states, we can 
always write 

f{x,Afp) = ln{n{xNp,Np))/K ■ (10) 

The function is nonzero only if < a; < 2. The symmetry 
([9]) implies that 

fix,Mp)^ f{2-x,J\fp) (11) 

In the statistical mechanics interpretation of the par- 
tition function (where /3 is an inverse temperature), 
f{x,J\fp) can be interpreted as a density of entropy. The 
existence of the infinite volume limit requires that 

limj^^^oof{x,JVp) = f{x) , (12) 

with f{x) volume independent. In the same limit, the 
integral ([5]) can be evaluated by the saddle point method. 
The maximization of the integrand requires 

f'{x) = (3 . (13) 

We believe that f{x) is strictly increasing for < x < 1 
with an absolute maximum at a; = 1. By symmetry, this 
would imply that f{x) is strictly decreasing for 1 < a; < 
2. We also believe that f'{x) is strictly decreasing and 
that Eq. (|13p has a unique solution (with positive /? if 
< a; < 1 and negative (3 ii 1 < x < 2). The numerical 
study of Sec. |TV] is in agreement with these statements, 
but we are not aware of mathematical proofs. Assuming 
that Eq. (|13|) has a unique solution, the infinite volume 
solution should be a; = P the average plaquette defined 
above. We can then convert an expansion for P into 
an expansion of /. If we want to include the volume 
dependence, the distribution has a finite width, and we 
should expand about the saddle point and perform the 
integration. In the following, we will work at large but 
finite volume and residual volume dependence in / will 
be kept implicit in equations. 

The behavior of f{x) for small x, can be probed by 
studying the model at large positive f3 (weak coupling 
expansion discussed in Sec. IVip . On the other hand, at 
small values of /3 (strong coupling expansion discussed 
in Sec. |V|, the partition function is dominated by the 
behavior of f{x) near its peak value a; = 1. For conve- 
nience, we introduce notations suitable for the study of 
the density of state near a; = 1 

g{y) = f{l + y). (14) 

g{y) is then an even function defined for — 1 < y < 1. 

IV. NUMERICAL CALCULATION OF n{S) 

To find n{S) numerically we will use a Monte Carlo 
simulation to create configurations of SU(2) for different 
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FIG. 2: Close-up of the patching process for 6*. 



values of /3. In the following example we will follow the 
steps we will use to find n{S) for a volume of 6^. We 
will start with 550 different sets of data ranging from 
P = 0.02 to /? = 11.00 in steps of 0.02 and with sizes of 
10^ configurations. To join the data from different values 
of f3 we will first create histograms of each set of data, 
each of these histograms is roughly Gaussian in shape. 
We then filter out the data that has statistics that are 
lower than half of the maximum bin. We can then remove 
the beta dependence by multiplying the height of each 
bin by e'^'^. We will be left with a series of arches which 
when overlayed on each other form the curve n{s). To 
create this overlay we will start with the lowest /3, which 
will correspond to the peak of n{s), and then take the 
logarithm of this. We will then look at the neighboring f3 
and do the same thing but then shifting it up or down so 
that the average distance in the bins overlapping with the 
first is zero. We will then continue in this manner until 
the supply of datasets has been exhausted. A portion 
of this process can be seen in Fig. O We then average 
the points for each bin together and divide both the bin 
width and height by Afp and shift the top of the curve to 
zero to make the final output, which can be seen in Fig. 
[3] for both 4^ and 6^. We see that they overlap well. 

We now consider the difference between two different 
volumes, as shown in Fig. 21 We can see in Fig. [5]that as 
we get closer to S/Afp — 1 this difference turns into noise, 
and as we get closer to S/Mp = we see a volume depen- 
dence growing. The results reported here correspond to 
the difference between 6'* and 4*. We have also studied 
the difference between 8** and 6"* and found consistent 
results. Calculation at larger volumes are much more 
computationally expensive and require many more sets 
of data because of the narrow width of the distributions. 
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SU(2) Patching Results 



SU(2) Diff 4* and 6* 





FIG. 3; Results of patching for 4'' and 6*. 



FIG. 6: The difference between ln{n{S))/J^p for 4* and 6* 
with f3 < 0.5 (close-up). 
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FIG. 4: The difference between ln{n{S))/Afp for 4^* and 6^. 
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FIG. 7: The difference between ln(n(5'))/7Vp for 4* and 6* 
divided by ln{S/Afp) 
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FIG. 5: The difference between ln{n{S))/JVp for 4"* and 6^ 
with /? > 0.5 (close up). 



V. STRONG COUPLING EXPANSION 

In this section, we discuss the strong coupling expan- 
sion of the logarithm of the density of state. We will work 
with the shifted function g{y) defined in Eq. The 
strong coupling expansion of P can be extracted from the 
expansion of InZ given in Ref. p^ . [isj using appropriate 
rescalings (for instance the /3 used there is one half of the 
(3 used here). The expansion is of the form 



P(/3) ^l+V] a2m-i/3' 



(15) 



m— 1 



The values of the coefficients are given in Table HI 

With periodic boundary conditions, the low order coef- 
ficients are volume independent. This can be understood 
from the exact translation invariance for the low order 
strong coupling graphs that provides a multiplicity that 
cancels exactly the 1 /Afp in Eq. ^ . Volume dependence 
may appear for graphs wrapping around the torus. The 
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TABLE I: Strong coupling expansion coefficients defined in 
tlie text. 
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simplest such graph is a straight line that closes into it- 
self due to the periodic boundary condition. It appears 
at order f3'^^ and has a reduced translation multiplicity 
since translation along the graph does not generate a 
new graph. This type of graphs produce 1/L corrections 
that to the best of our knowledge have not been studied 
quantitatively. In the following, we will ignore such ef- 
fects, but a study of the contribution of graphs with a 
nontrivial topology would certainly be interesting. 
We will plug the expansion of P in the expansion 

g{y) ^ 32™y"" ■ (16) 

m=0 

At lowest order we have y ~ ai/3 and the saddle point 
Eq. yields 2g2y — 2g2aif3 ~ /3 which implies 52 = 

l/(2ai). This procedure can be followed order by order 
in p. The results are shown in Table HI 

Since n{S) is zero for S* = and 2Afp, we expect loga- 
rithmic singularities at x = and 2 for f{x) and y = ±1 
for g{y). This singularities will cause the strong coupling 
series to diverge when \y\ > 1. Consequently, we define 
the subtracted function 

hiy)^giy)^Ailnil^y')) . (17) 

The coefficient A will be calculated using the weak cou- 
pling expansion in Sec. IVII In the infinite volume limit, 
we have A — 3/4. Expanding 

hiy) c J2 h2my'"' , (18) 

we obtained coefficients that are showed in Table U for 
A = 3/4. The coefficients g2m and ft,2m are also shown on 
a logarithmic scale in Fig. [H This graph shows that the 
two types of coefficients become rapidly of the same or- 
der, which indicates singularities in the complex y plane 
for smaller values of \y\ than the ones at ±1. 

In Fig. [9l we show the error made at successive order 
of the strong coupling expansion of the plaquette. We 
then show successive approximation of f{x) (Fig. 
and the corresponding errors (Fig. Ilip . 



FIG. 8: Logarithm of the absolute value of g2m and /12 
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FIG. 9: Logarithm of the absolute value of the difference be- 
tween the numerical data and the strong coupling expansion 
of P at successive orders. For reference, we give the estimated 
numerical error on P. 
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FIG. 10: Numerical value of f[x) compared to the strong 
coupling expansion at successive orders. 
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SU(2) 6 strong coupling 
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FIG. 11: Logarithm of the absolute value of the difference be- 
tween the numerical data and the strong coupling expansion 
of / at successive orders. For reference, we give the estimated 
numerical error on /. 

We can now compare the apparent convergence of P 
and /. From Fig. [S] we see that the larger order errors 
cross between p = 1.5 and 2. For values of /3 larger, in- 
cresing the order increases the error. This is the sign of a 
finite radius of convergence [l3| • Similarly, the larger or- 
der errors for / cross for x between 0.5 and 0.6 which are 
approximately the values of P in the (3 interval of cross- 
ing. Consequently, it seems like the convergence proper- 
ties of the two expansions are the same (finite radius of 
convergence) . 

VI. WEAK COUPLING EXPANSION 

In this section, we discuss the weak coupling expansion 
of f{x). The starting point is the expansion of P in 
inverse powers of /3 

m— 1 

We then assume the behavior 

fix) ~ A ln{x) + Z™^" ■ (20) 

m=0 

Using the saddle point Eq. (fT5|) . and using the leading 
large (3 and small x terms we find 

/? ~ A/a; ~ A/(6i//3) , (21) 

which implies A = bi dX infinite volume. The procedure 
can be pursued order by order without difficulty. The 
result for the two lowest orders is 

/i = h^/bi 

f2 = ib3b,-bl)/{2bl) 



Numerical experiments indicate that the two series have 
the same type of growth (power or factorial) . Note that 
/o cannot be fixed by the saddle point equation. The 
overall height of / depends on the behavior near x — 1 
(if we insist on normalizing n{S) as probability density) 
and it seems unlikely that it can be found by a weak 
coupling expansion. 

At finite volume, the saddle point calculation of P 
should be corrected in order to include 1/V effects {V 
the number of sites , for a symmetric lattice). If we 
perform the Gaussian integration of the quadratic fluc- 
tuations, and use the V dependent value of bi given in 
Eq. 123]) below, we find after a short calculation that the 
coefficient A of ln{x) is 

A = (3/4)-(5/12)(l/F) . (22) 

This leading coefficient correction, predicts a difference 
of — 0.0013/n(x) for the difference between f{x) for a 4** 
and 6^ and is roughly consistent with Fig. [71 

Our next task is to find the values of b„i. A closed 
form expression can be found [H, [IB] for bi. For the case 
Nc = 2 and _D = 4, we obtain 

6i = (3/4)(l- 1/(31-)) . (23) 

The (— 1/(3F)) comes from the absence of zero mode 
{—1/V) in a sum calculated in [l^ plus the contribu- 
tion of the zero mode with periodic boundary conditions 
(-1-2/(31^)) calculated in [l^|. Numerical values for 62 can 
be found in Ref. [lB| and for 63 in Ref. [l3| • In these 
Refs., several sums are calculated numerically at partic- 
ular volumes that do not include 6*. Rough extrapo- 
lations from the existng data indicate that for V — G'^ 
uncertainties are less than 0.0002 for 62 and 0.0008 for 
63. For /3 > 3, these effects are close to the numerical 
errors for P. In the following, we use the approximate 
values 62 = 0.1511 and 63 = 0.1427 for V = 6^. 

We are not aware of any calculation of bm for m > 4 
for SU(2) In the case of SU{3), calculations up to or- 
der 10 [igl and 16 [l^ are available and show remarkable 
regularities. Using the assumption Q that dP/dfi has a 
logarithmic singularity in the complex (3 plane and inte- 
grating, we obtained [20| the approximate form 

Y bmr" - C(Li2(/3-V(/3™' + *r)) + h.c , (24) 

m— 1 

with 

Li2(a;) =^xVfc^ . (25) 

k=l 

We believe that at zero temperature, the new parameter 
r which measures the (small) distance from the singu- 
larity to the real axis in the 1//3 plane stabilizes at a 
nonzero value in the infinite volume limit. For reasons 
not fully understood, this parametrization of the series 
turns out to work very well for SU{3). For instance, by 
fixing the value of F in the middle of the allowed range 
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and using the values of 5g and 610, we obtain values of 
the lower order coefficients with a relative error of 0.2 
percent for &§ a-nd that increases up to 5 percent for 63. 
In the limit F = 0, the parametrization provides simple 
predictions for instance 63/62 — (4/9)/3m- The location 
of the Fisher's zeros for S'C/(2) [| suggests /3,„ = 2.18. 
This implies 63/62 — 0.969 in good agreement with our 
numerical estimate 63/62 ~ 0.944. In the following we 
use the values = 2.18, F = 0.18//3^ ~ 0.038 (see [|) 
and we fixed C = 0.0062 in order to reproduce 63. The 
numerical values of 6^ and the corresponding values of 
fm are displayed in Table HIl 

We have compared the weak coupling expansion of P 
with numerical values in the case V = 6'^. The results 
are shown in Fig. 1121 In the region where the curves are 
smooth, the error decrease with the order and appears to 
accumulate. This is very similar to the case of SU (3)[20l|. 
However, it is clear that more reliable estimates for m > 4 
would be desirable for SU (2). It should be noted that for 
large f3, the noise in the error is at the same level as the 
numerical error on P. This would not be the case if we 
had not included the contribution of the zero mode to 61 
as shown in the second part of Fig. [121 

We have compared the weak coupling expansion of 
f{x) with numerical values in the case = 6^. The re- 
sults are shown in Fig. [131 The differences are resolved 
in Fig. [131 In these graphs we have taken /o = —0.14663 
which maximizes the length of the accumulation line on 
the left of Fig. [H 
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0.3643 


10 


2.9652 


0.5883 


11 
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14 


33.690 


5.3207 


15 


63.702 


9.6945 



TABLE II: Weak coupling coefficients defined in Sec. |Vll The 
choice of fei corresponds to = 6'' . 
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FIG. 12: Logarithm of the absolute value of the difference 
between the numerical data and the weak coupling expansion 
of P at successive orders (above). For reference, we give the 
estimated numerical error on P. The graph below is the same 
except that we have not included the zero mode in bi. 
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FIG. 13: Numerical value of f{x) compared to the weak 
coupling expansion at successive orders. 
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numerical error on /. 



VII. EXPANSION IN LEGENDRE 
POLYNOMIALS 

We now consider the function h(y), which is g{y) with 
the logarithmic singularity subtracted as defined in Eq. 
P?)l . This is a bell shaped even function defined on the 
interval [—1,1] and shown on Fig. We can expand 
this function in terms of the even Legendre polynomials. 

Hy) = 1^rr.P2miy) (26) 

m=0 

The 92m can be determined from the orthogonality re- 
lations with interpolated values of h{y) to perform the 
integral. A minor technical difficulty is that we do not 
have numerical data all the way down to y = — 1. This is 
because as y ^ — 1+, or in other words a; — > 0"*", /? — > +00 
where the plaquette distribution becomes infinitely nar- 
row. Consequently there is a small gap in the numerical 
data that needs to be filled. Fortunately, this is precisely 
where the weak coupling expansion works well. Using 
the weak coupling expansion (including the overall con- 
stant), subtracting Aln(x{2 — x)) and shifting to the y 
coordinate, we obtained the approximate behavior near 
y = — 1 for the B** data: 

h{y) ~ 0.2145 + 1.2961?; + 0.5261^^ + 0.1109^^ (27) 

In order to estimate the error associated with this ap- 
proximation we have compared with an extrapolation of 
a quadratic fit of the leftmost part of the data. In order 
to give an idea of the volume effects, we have also used 
the second method on a 4'* lattice. The results are shown 
in Table Hill This indicates that the variations are small, 
increase with the order in relative magnitude and that 
the volume effects are stronger than the dependence on 
the extrapolation procedure. The logarithm of the coeffi- 
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FIG. 15: Legendre polynomial coefficients q2m with the three 
methods described in the text. 

cients is shown in Fig. [TSlwhich illustrate the exponential 
decay of the coefficients. 

The expansion provides excellent approximation of 
h(y) shown in Fig. 1161 The errors are resolved in Fig. 
[T71 It is also possible to calculate P{f3) by solving the 
saddle point Eq. (|13p using successive approximations 
for h. This is shown in Figs. [18] and [191 The spikes 
in the error graphs correspond to change of sign of the 
errors. It is important to notice that the quality of the 
approximations improves with the order in all region of 
the interval. 



VIII. CONCLUSIONS 

We have calculated the density of states for SU{2) lat- 
tice gauge theory. The intermediate orders in weak and 
strong coupling agree well in an overlapping region of 
action values as shown in Fig. [201 However, the large or- 



method 


4" +m 


6'^ +&t 


6^+ ^ 


m 




92 m 







-0.30034 


-0.30095 


-0.30096 


1 


-0.47963 


-0.48159 


-0.48164 


2 


0.1488 


0.14853 


0.14845 


3 


-0.03215 


-0.0309 


-0.03099 


4 


-0.00822 


-0.00843 


-0.00852 


5 


0.01156 


0.01114 


0.01107 


6 


-0.00363 


-0.00305 


-0.00308 


7 


-0.00186 


-0.00179 


-0.0018 


8 


0.00194 


0.00146 


0.00147 


9 


0.00008 


0.00026 


0.00028 


10 


-0.00094 


-0.00069 


-0.00067 



TABLE III: Legendre polynomial coefficients q2m with the 
three methods described in the text. 
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FIG. 16: h{y) together with the expansion in Legendre poly- 
nomials up to order 20. 



FIG. 19: Logarithm of the absolute value of the difference 
between the numerical data for P and expansions in Legendre 
polynomials at successive orders 
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FIG. 17: Logarithm of the absolute value of the difference be- 
tween the numerical data for h{y) and expansions in Legendre 
polynomials at successive orders 
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FIG. 18: P together with the expansion in Legendre polyno- 
mials up to order 20. 



SU(2) weak and strong coupling 
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FIG. 20: Weak and strong coupling expansion of / at a few 
intermediate orders. 



der behaviors of these expansions appear to be similar to 
the corresponding ones for the plaquette. Volume effects 
can be resolved well for small actions values. Correc- 
tions to the saddle point estimate need to be developed 
systematically. Aprroximation of a subtracted quantity 
by Legendre polynomials looks very promising and works 
well uniformly. We plan to use this approximate form to 
look for Fisher's zeros. 

The density of states can be calculated in more general 
situations. For instance, 



/■2A/'p 

ZiP,m)^ d5n(5,{A})e-''^ , (28) 

Jo 
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with 



n{S,{(3^}) = (29) 
YlJdUiSiS - ^(1 - (l/Ar)i?err([/p))) (30) 



'E./3.(i-x.(;7p)M) 



(31) 



and Xi ^ complete set of SU{2) characters. This is a 
type of action which naturally arises in the RG studies 
of SU{N) lattice gauge theories. It is possible to ap- 
ply exact renormalization group transformation [T, '5| or 
the MCRG procedure ]21] to the partition function in 
order to define the couplings. Following the analogy be- 
tween f — (3 and V = J for the effective potential V 



in presence of a source J in scalar models, it would be 
interesting to study finite size effects from this point of 
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